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ERROR ESTIMATES AND A TWO GRID SCHEME FOR 
APPROXIMATING TRANSMISSION EIGENVALUES* 
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Abstract. In this paper, using the linearization technique we write the Helmholtz transmission 
eigenvalue problem as an equivalent nonselfadjoint linear eigenvalue problem whose left-hand side 
term is a selfadjoint, continuous and coercive sesquilinear form. To solve the resulting nonselfadjoint 
eigenvalue problem, we give an conforming finite element discretization and establish a two grid 
discretization scheme. We present a complete error analysis for both discretization schemes, and 
theoretical analysis and numerical experiments show that the methods presented in this paper can 
efficiently compute real and complex transmission eigenvalues. 
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1. Introduction. Transmission eigenvalues not only have wide physical applica¬ 
tions, for example, they can be used to obtain estimates for the material properties of 
the scattering object HE], but also have theoretical importance in the uniqueness and 
reconstruction in inverse scattering theory m- Many papers such as [a [131 m [Ml [29] 
study the existence of transmission eigenvalues, and [211115] explore upper and lower 
bounds for the index of refraction n{x) from knowledge of the transmission eigenval¬ 
ues. The attention from the computational mathematics community is also increasing 

(see, e.g., 11 [HillEIllMl 1131211132 )■ 

In this paper, we consider the Helmholtz transmission eigenvalue problem; Find 
k € C, w,a £ L^{D), w — a £ H^{D) such that 


(1.1) 

Aw + k^n{x)w = 0, 

(1.2) 

Act + fc^CT = 0, in 

(1.3) 

w — a = 0, on dD 


dw da 

(1.4) 

dv dv 


where D {d = 2, 3) is a bounded simply connected set containing an inhomoge¬ 

neous medium, and v is the unit outward normal to dD. 

From im we know that for u = w — a £ Hq{D), the weak formulation for the 
transmission eigenvalue problem (HI)-(El can be stated as follows: Find £ C, 
k^ ^ 0, u £ IIq(£>)\{0} such that 


(1.5) 


n(x) — 1 


(Au + k^u), Av + k'^n{x)v)o = 0, 




where (V', </5)o = Jo4’^dx denotes the L‘^{D) inner product. We denote A = as 
usual, then (|1.5I1 is a quadratic eigenvalue problem. 

In recent years, based on this weak formulation, various numerical methods have 
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been presented to solve the transmission eigenvalue problem. The first numeri¬ 
cal treatment appeared in |14j where three finite element methods were proposed 
for the Helmholtz transmission eigenvalues which have been further developed in 
[B [9l [HI im |22l El] • Among them [Ml [Hi El] studied the conforming finite ele¬ 
ment method, iniDjiiii] mixed finite element methods, [1] spectral-element method 
and [19j the Galerkin-type numerical method. Inspired by these works, this paper 
further studies the conforming hnite element method and has three features as 
follows: 

(1) A complete error analysis is presented. Due to the fact that the problem is 
neither elliptic nor self-adjoint, once its error analysis was viewed as a difficult task. 
Sun [31] first proved an error estimate for the conforming finite element approxi¬ 
mation of m- Following him Ji et al. [22] presented an accurate error estimate, but 
Theorems 1-2 therein are only valid for real eigenvalues and rely on the conditions 
of Lemma 3.2 in m which are not easy to verify in general, especially for multi¬ 
ple eigenvalues. In this paper, we use the linearization technique in [33] to write the 
weak formulation (lES]) as an equivalent nonselfadjoint linear eigenvalue problem (|2.6I1 . 
Then the conforming finite element approximation eigenpair {Xh,Uh) of (12.611 is 
exactly the one of (lES]). Fortunately, the left-hand side term of (12.611 is a selfadjoint, 
continuous and coercive sesquilinear form, thus we can use Babuska-Osborn’s spectral 
approximation theory |2] to give a complete error analysis for the conforming finite 
element approximation of (12.611 . namely the error analysis of (11.511 . The difficulty of 
the error estimates in low norms lies in the nonsymmetry of the right-hand side term 
of (12.611 that involves derivatives. We introduce two auxiliary problems and overcome 
this difficulty by the Nitche technique in a subtle way. Our theoretical results are 
proved under general conditions and are valid for arbitrary real and complex eigen¬ 
values. 

(2) A finite element discretization with good algebraic properties is presented. 
We give an conforming finite element discretization (see (13.111 or (15.111 1 that will 
lead to a positive definite Hermitian and block diagonal stiff matrix. Then we use 
this discretization to solve the transmission eigenvalue problem numerically and ob¬ 
tain both real and complex transmission eigenvalues of high accuracy as expected. 
Similar discretizations already exist in the literatures. Colton et al. [M] established 
the conforming finite element discretization for the formulation (II.5p (see (4.3) in 
[14] 1. by starting with discretizing (IE3 into a quadratic algebra eigenvalue problem 
and then linearizing it, which is in the opposite order of our treatment. Though our 
discrete form m is formally different from the one in [14] for the resulting mass ma¬ 
trix in [14] is positive definite Hermitian and block diagonal, the two discrete forms 
are equivalent from the linear algebra point of view. Gintides et al. m has used 
square root of matrices, which is a standard tool, and a Hilbert basis in Hq{D) as 
test functions to discretizite (11.51) into a linear eigenvalue problem (see (2.14) in [19] 1 
and proved the convergence for the approximate eigenvalues. 

(3) A two grid discretization scheme is proposed. As we know, it is difficult to solve 
nonselfadjoint eigenvalue problems in general. So, Sun m adopted iterative methods 
to work on a series of generalized Hermitian problems and finally to solve for the real 
transmission eigenvalues. Ji et al. [22] combined the iterative methods in [31] with 
the extended finite element method to establish a multigrid algorithm for solving the 
real transmission eigenvalues. This paper establishes a two grid discretization scheme 
directly for the nonselfadjoint eigenvalue problem (12.611 which is suitable for arbitrary 
real and complex eigenvalues. The two grid discretization scheme is an efficient ap- 
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proach which was first introduced by Xu |34j for nonsymmetric or indefinite problems, 
and successfully applied to eigenvalue problems later (see, e.g., [T71 [Ml [3S113S113Z] and 
the references therein). With the two grid discretization scheme, the solution of the 
transmission eigenvalue problem on a fine grid nh is reduced to the solution of the 
primal and dual eigenvalue problem on a much coarser grid -kh and the solutions of 
two linear algebraic systems with the same positive definite Hermitian and block di¬ 
agonal coefficient matrix on the fine grid tt^i, and the resulting solution still maintains 
an asymptotically optimal accuracy. The key to our theoretical analysis is to find a 
{u*u,u}*u) and prove that \B[{uht^jJh)i{u*u,uj*u))\ (see Step 1 of Scheme 4.1) has a 
positive lower bound uniformly with respect to the mesh size H, which is also critical 
for further establishing multigrid method and adaptive algorithm. We successfully 
apply spectral approximation theory to solve this problem (see Lemma 4.1 and Re¬ 
mark 4.1). 

The conforming finite element discretization is easy to implement under the 
package of iFEM m with MATLAB, and the numerical results indicate that this dis¬ 
cretization is efficient for computing transmission eigenvalues. Moreover, we can fur¬ 
ther improve the computational efficiency by using the two grid discretization scheme. 
Regarding the basic theory of finite element methods, we refer to [am [121127]. 
Throughout this paper, C denotes a positive constant independent of the mesh 
size h, which may not be the same constant in different places. For simplicity, we use 
the symbol a "^b to mean that a < Cb. 

2. A linear formulation of (|1.5I) . In this paper, we suppose that the index of 
refraction n S L°°{D) satisfies the following assumption 

1 -I- 5 < inf n(a;) < n(x) < supn(a;) < oo, 

D D 

for some constant 5 > 0, although, with obvious changes, the theoretical analysis in 
this paper also holds for n strictly less than 1. 

Inspired by the works in [9i[ni[2T], we now apply the linearization technique in 
[33j to write the weak formulation (11.51) as an equivalent linear formulation. 

From dUD we derive that 

(—^Au, Av)o + k'^{ —Au)o 

77 ,- 1 77,-1 

(2.1) -|-fc^(A-», ^ v)o + k'^{ ^ u,v)o = 0, Vue Hq{D). 

n — 1 n — 1 

Introduce an auxiliary variable 

( 2 . 2 ) uj = k^u, 

then 

(2.3) (w, z)o = fc^(u, z)o, Vz e L^{D). 

Thus, combining (12.11) and (12.3p . we arrive at a linear formulation: Find e C and 
nontrivial {u,uj) € Hq{D) x L‘^{D) such that 

( ^ T Au,Av)o = -k'^{—^u,Av)o 

77 ,- 1 77 —1 

77 77 

(2.4) -k^{Au, - vu)o-fc2( -ra;,u)o, Vue 

77 ,- 1 77,-1 

(2.5) (w, z)o = k'^{u, z)o, Vz e L'^{D). 
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This is a nonselfadjoint linear eigenvalue problem. 

Suppose that {k'^,u) is an eigenpair of (II.5L then it is obvious that is 

an eigenpair of (I2.4D - (I2.5I) . On the other hand, suppose that satisfies (12.4^ - 

(12.51) . from (12.51) we get w = and substituting it into (12.41) we get (l2.1l) for m)- 
The above argument indicates that (j2.4l) - (l2.5l) and (11.51) are equivalent. 

Let H~^{D) be the “negative space”, with norm given by 


ll/ll-i= sup 


Ml 


1 = 1 , 2 . 


Define the Hilbert space H = Hq{D)xL^{D) with norm ||(u,a;)||H = l|u ||2 +||w||o: 
and define H'’ = H'^{D) x with norm ||(u,w)||h‘» = ||m||s+ ||w||s-2, s = 0,1, 

= H. It’s obvious that H ^ H'* (s = 0,1) compactly (see[5l[35] ). 

Denote A = Let 

A{{u, w), {v, z)) = (^—Am, Ai;)o + (w, z)o, 
n — 1 

1 Ti 71 

B{{u,uj), iv,z)) = -(--u, A?;)o - (Ait,--u)o - (-rW,u)o + (u,z)o, 

n — 1 n — 1 77,-1 

then (I2.4l) - p.5l) can be rewritten as: Find A G C, (u,a;) G H\{0} such that 


( 2 . 6 ) 


A{{u,uj), {v,z)) = XB{{u,u}), iv,z)), \f{v,z) G H. 


Thus we get the following theorem. 

Theorem 2.1. If {k^,u) is an eigenpair of and w = k'^u, then (kf,u,uj) 

satisfies and if {k'^ ,u,uj) satisfies 12.61) . then {k'^,u) is an eigenpair of 11.5]) and 

uj = k'^u. By calculation we derive for any (u,a;), {v, z) G H, 


1 


1 


A((u,uj), (v, z)) = / - AuAv + uizdx= / - AvAu + zuidx = A((v, z), (u,ui)), 

J n-l J n-1 

D D 

\A{{u,u),{v,z))\ < ||Au||o||Au||o + ||w||olkllo < ||(u, w)||h||(u,2:)||h, 

1 


A{{u,uj), iu,uj)) = 


n — 1 


AuAu + ujijjdx > ||(u,w)"^ 




D 


i.e., A{-, •) is a selfadjoint, continuous and coercive sesquilinear form on H x H. 

We use A{-, •) and || • m = A{-, ■)^ as an inner product and norm on H, respec¬ 
tively. Then || • m is equivalent to the norm || • ||h- 
When n G L°°{D), \/{f,g), {v,z) G H, we deduce 

1 71 71 

\B{{f, g), {v, z))\ = \- (--/, A?;)o - (A/,--u)o - (- -g, v)o + (/, z)o| 

77 — 1 77—1 77—1 

1 77 77 

< II-7/llolkl|2 + II--A/||_i||u||i + II --5||-l||u||i + ||/||o||z||o 

77—1 77—1 77—1 

(2.7) < (ll/llo + II^A/||_i + ||^g||_0||(u,z)||H. 

77 — 1 77 —1 

When n G y{f,g) G H^, V(z;,z) G H, we have 

1 77 77 

\B{{f, g), {v, z))| = |(V(--/), V^;)o + (V/, V(- -v))o - (- -g, v)o + if, ^)o| 

77—1 77—1 77—1 

<ll/llilklli + ll/llilklli + llffll-ilklli + ll/llill^ll 

(2.8) <||(/,5)IIhi||(^^,^)||h. 


-1 
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And when n G y{f,g) G H°, y{v,z) G H, we have 

1 71 T) 

\B{{f,g), {v,z))\ = I - (- -f,Av)o - (/, A(- -v))o - (- -g,v)o + {f,z)o\ 

n — 1 n — 1 n — 1 

< Il/ll0|k||2 + ||/||0|k||2 + ||g||-2|k||2 + ||/||0||^||0 

(2.9) < ||(/,ff)||Ho||(l^,^)||H. 

We can see from (|2.7I1 - (I2.9D that for any given {f,g) G H'* (s = 0,1, 2), B{{f,g), {v, z)) 
is a continuous linear form on H. 

The source problem associated with (12.61) is given by: Find {ip, </?) G H such that 

(2.10) A{{ip,(p),{v,z)) = B{{f,g),{v,z)), V(u,z)gH. 

From the Lax-Milgram theorem we know that (12.101) admits a unique solution. There¬ 
fore, we define the corresponding solution operators T : H® —H (s=0,l,2) by 

(2.11) A{T{f,g),{v,z)) = B{{f,g),{v,z)), V(u, z) G H. 

Then (12.61) has the equivalent operator form: 

(2.12) T{u,uj) = \-^{u,<jj). 

Namely, if {X,u,u}) is an eigenpair of (12.6L then {X,u,uj) is an eigenpair of (j2.12l) . 
Conversely, if {X,u,uj) is an eigenpair of (I2.12p . then {X,u,uj) is an eigenpair of (12.61) . 

Theorem 2.2. Suppose n G L°°{D), then T : H —>■ H is compact; suppose 
n G 1F^“®’°“(I1) (s = 0,1), then T : H® —>■ H® is compact. 

Proof. When n G L°°{D), let (u, z) = T{f,g) in (I2.11L then from (12.71) we have 

lin/,ff)ll^ < A{T{f,g),T{f,g)) = B{{f,g),T{f,g)) 

71 71 

< (ll/llo + II^A/II-I + ||;^5ll-i)lin/,ff)llH, 

thus 

(2.13) ||r(/,g)||H< ll/llo+ ||^A/||_i + ||^5||_i, V(/,g)GH. 

n — 1 n — 1 

Let S = X ^2 be any bounded set in H, because of the compact embedding 
Hq{D) ^ L‘^{D) and L‘^{D) ^ H~^{D), then is relatively compact in Lf{D), 

{v : V = .^^A/, / G S'!} is relatively compact in H~^{D) and {v : v = G «S' 2 } 

is relatively compact in H~^{D). And from (I2.13|) we conclude that T : H —>■ H is 
compact. 

When n G VF^“®’°°(Z1) (s = 0,1), in (12.111) . let {v,z) = T{f,g), then from (12.81) 
and (12.91) we have 

mLgWu < A{T{f,g),T{f,g)) = B{{f,g),T{f,g)) < \\{f, g)\\ii4nf, g)\\H, 

thus 

(2.14) ||T(/,5)||h< ||(/,5)||h», V(/,5)gH® (s = 0,1), 

which implies that T : H® H is continuous. Due to the compact embedding 
H ^ H®, T : H® —H® is compact. □ 

Consider the dual problem of (12.61) : Find A* G C, {u*,u)*) G H\{0} such that 

(2.15) A((u,z),(w%a;*))=FB((u,z),(u*,cc*)), V(i;,z)gH. 
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Define the corresponding solution operators T* : H® —^ Hby 

(2.16) A{{v,z),T*{f,g))=Bi{v,z),{f,g)), V(z;,z)gH. 

Then (12.151) has the equivalent operator form: 

(2.17) T*{u*,io*) = X*-\u*,uj*). 

It can be proved that T* is the adjoint operator of T in the sense of inner product 
A{-, ■). In fact, from (j2.11|) and ()2.16ll we have 

{v,z)) = B{{f,g), {v,z)) 

(2.18) =Ai{f,g),T*{v,z)), V(/, 5 ), (u, z) G H. 

Hence the primal and dual eigenvalues are connected via A = A*. 

3. The conforming finite element approximation and its error esti¬ 
mates. Let TTh be a shape-regular grid of D with mesh size h and C H^{D) be a 
piecewise polynomial space defined on tt/j ; for example, is the finite element space 
associated with one of the Argyris element, the Bell element or the Bogner-Fox-Schmit 
element (BFS element) (see [H]). Thanks to (12.21) . we choose H^, = x S^. Then 
H;, C H be a conforming finite element space. For the three finite element spaces 
mentioned above, since U?i>o dense in both L‘^{D) and Hq{D), it is obvious 

that the following condition (Cl) holds: 

(Cl) If ('0, (/?) G H, then as h —>■ 0, 

inf !!(?/',(/?)-(u, 2 ;)IIh ^ 0 . 

(t),z)6Hh 

From the operator interpolation theory (see e.g. i), the following (C2) is also valid: 
(C2) If ^ G H^{D) n H'^+^{D) (r G (0, 2]), then 

inf IIV’- u||s < s = 0 , 1 , 2 . 

«eS'* 

Throughout this paper, we suppose that (Cl) is valid. 

The conforming finite element approximation of (12. 6 p is given by the following: 
Find Xh G C, {uh,ujh) G H^\{0} such that 

(3.1) A{{uh,ujh),{v,z)) = XhB{{uh,iZ}h),{v,z)), V(u,z)GH?i. 

Consider the approximate source problem: Find G H/j\{0} such that 

(3.2) A{{iph,‘Ph),{,v,z)) = B{{f,g),{v,z)), V(-i;,z)gH,,. 

We introduce the corresponding solution operator: Th : H® —H/j (s=0,l): 

(3.3) A{n{f,g),{v,z)) = B{{f,g),{v,z)), V(u,z)gH;,. 

Then (13.11) has the operator form: 

(3.4) T)i(u/j, cc/j) = A^ ^h)• 

Define the projection operators : Hq{D) —)• and : L‘^{D) —>• by 

(^--A(u-Au)o = 0, VuGS”^, 
n — 1 

(w — P^u}, z)o = 0, Vz G S^. 


(3.5) 

(3.6) 
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Let 

Ph{u,uj) = {Plu,Pluj), V(M,w)eH. 

Then : H —H/j, and 

A((m,w) - Phiu,u;), {v,z)) = A{{u,u;) - (P^u^P^uj), {v,z)) 

(3.7) = {^—-A{u- pIu),Av)o +{uj - pIuj,z)o = 0, \/{v,z)€lih, 

n — I 

i.e., is the Ritz projection. 

Combining (13.71) . (12.111) with (13.31) . we deduce for any (u,w) G H that 

A{PhT{u, to) - Th{u, oj), (v, z)) = A(PhT(u, w) - T(u, oj), (v, z)) 

+ A(T(u, uj) - Th{u, uj), {v, z)) = 0, V(n, z) G H/^, 

thus we get 

(3.8) n = PhT. 

Theorem 3.1. Let n G L°°{D), then 

(3.9) ||T-T^,||h->0, 
and let n G W^~^’°°{D) , then 

(3.10) ||T-r^||H» ^0, s = 0,l. 


Proof. When n G L°°{D)^ for any (/, 5 ) G H, from (Cl) we have 
II(/ - P,)T(/,g)||H = ||(T(/,g) - P,r(/,ff)||H < ^ inf \mf,g) - (u,2)||h -> 0. 

Since T : H — 7 > H is compact, T{{f,g) G H, ||(/,g)||H=i} is relatively compact. Thus 
by the definition of operator norm we have 

r-T,||H= sup ||(T-r^){/,g)||H 

(/.s)eH,||(/,g)||H = l 

sup ||(/-P^)r(/,5)||H^0. 

(/,g)6H,||(/.g)||H = l 

When n G W^“®’°“(Z1), for s = 0,1, we have 

||(/-p^)r(/,5)||H= < ||(/-p^)r(/,5)l|H ^ 0 . 

Since T : H® ^ H® is compact, T{{f,g) G H®, ||(/,5)||h'»=i} is relatively compact, 
thus we have 

||r-T,||H«= sup ||(T-r^)(/,g)||H» 

(/.g)eHM|(/.g)||H«=l 

sup ||(/-P^)r(/,5)||H» ^0. 

(/.g)eHM|(/.g)||H<.=l 


The proof is completed. □ 
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The conforming finite element approximation of (j2.15l) is given by: Find G C, 
K,cc,t)GHA{ 0 } such that 

(3.11) = V(n,z)GH^. 

Define the solution operator : H® —>• Hfi satisfying 

(3.12) A{{v,z),mf,g)) = B{{v,zl{f,g)), y{v,z)GUn. 

(13.111) has the following equivalent operator form 

(3.13) niulu:l) = Xr\ulc,l). 

It can be proved that is the adjoint operator of Th in the sense of inner product 
A{-, •). In fact, from (13.31) and (13.121) we have 

MTh{f,g), (v^z)) = B{{f,g), {v,z)) 

(3.14) =Aiif,g),mv,z)), V(/,g),(u,z)GH^. 

Hence, the primal and dual eigenvalues are connected via Xh = A)). 

We need the following lemma (see Lemma 5 on page 1091 of [TH]). 

Lemma 3.2. Let HT/j — T||h —0. Let {Aj} be an enumeration of the eigenvalues 
ofT, each repeated according to its multiplicity. Then there exist enumerations {Aj^/j} 
of the eigenvalues ofTh, with repetitions according to multiplicity, such that Xj^h 

A, (j > i;. 

In this paper we suppose that {Ajj and {Aj_ft,} satisfy the above lemma, and let 
A = Afe be the kth eigenvalue with the algebraic multiplicity q and the ascent a, 
Afc = Afc+i = • • •, Xk+q-i- Since ||T,i - THh ^ 0, q eigenvalues Xk,h, ■■■, Xk+q-i,h of 
will converge to A. 

Let E be the spectral projection associated with T and A, then R{E) = N{{X~^ — 
T)“) is the space of generalized eigenvectors associated with A and T, where R denotes 
the range and N denotes the null space. Let Eu be the spectral projection associated 
with Th and the eigenvalues Xk,h, ■ • ■, Xk+q-i,h', let Mh{Xi^h) be the space of gener¬ 
alized eigenvectors associated with Xi^h and Th, and let Mh{X) = J2i=k~^ Mh{Xi^h), 
then R{Eh) = Mh{X) if h is small enough. In view of the dual problem ()2.15l) and 

(13.111) . the definitions of E*, R{E*), Ef^, Mh{X*ff), Mh{X*) and R{Eff) are analogous 
to E, R{E), Eh, Mh{Xi^h), Mh{X) and R{Eh) (see [2]). 

Given two closed subspaces R and U, denote 


5{R,U) = 

II 

e{R, u)s = 


sup iiif ||(■u,a;) - (w,2)||h, 

{v,z)&U 

(«,-)IIh = i 

sup inf ||(u,a;) - (■i;,2;)||H-, 

{v,z)GU 

II(«,<*')IIhs = i 


s = 0,1. 


We define the gaps between R{E) and R{Eh) in || . ||h as 

S{R{E),R{Eh)) = max{J(i?(F;), RiEh)), S{R{Eh), RiE))}, 


and in || . ||h« as 

9iRiE),R{Eh)), = ma^{9iRiE),R{Eh))s,0{R{Eh),RiE))s}. 
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Define 


eh{X) = sup inf ||(u,a;) - (u, 2 :)||h, 

(u,u)eii(E) 

1I(u,^)|Ih=i 


4(A*) 


sup inf ||(u*,w*) - (u,z)||h- 

(u*,u*)eH(E*) («, 2 )eHfc 


It follows directly from (Cl) that 

e?,(A) ^ 0 (h ^ 0), 4(A*) 0 (h ^ 0). 

Suppose that R{E),R{E*) C Hn (r G (0,2]), then from (C2) we get 


(3.15) e^(A) < el(A*) < hC 

Further suppose that R{E),R{E*) C (i?®(D))^, C Hq{D) is the Argyris finite 
element space, then from the interpolation theory we have 

(3.16) eh{X) < h\ eliX*) < h\ 

Note that when the functions in R{E) and R{E*) are piecewise smooth (13.151) and 
(13.161) are also valid. 

Thanks to [2], we get the following Theorem 3.3. 

Theorem 3.3. Suppose n G L°°{D), then 

(3.17) S{R{E),RiEh))<eh{X), 

k+q-l 

(3.18) |A-(- ^ A-i)-i|<e,(A)4(A*), 

® j=k 

(3.19) |A — Xj^h\ ^ [e;i(A)e)](A*)] “ , j = k, k + 1, ■ ■ ■, k + q — 1. 

Suppose {uh,u!h) with ||(M^,w^i)m = 1 is an eigenfunction corresponding to Xj^h (j = 
k,k + 1^ ■ ■ ■ ,k + q — i), then there exists an eigenfunction (u, ui) corresponding to X, 
such that 


(3.20) 


\\{uh,ujh) - (u,a;)||H < e/i(A)“. 


Proof. From Theorem 3.1 we know IjT — T;,||h ^"0 (h^-O), thus from Theorem 
8.1, Theorem 8.2, Theorem 8.3 and Theorem 8.4 of [2] we get the desired results 
(I3.17p . (I3.18|) . (13.191) and (I3.20|) . respectively. □ 

Next we discuss the error estimates in the || • ||h* (s = 0,1) norm by using the 
Aubin-Nitsche technique. 

We need the following regularity assumption: 

R(D). For any ^ G H ®(iA) (s = 0,1), there exists if G 77^+’'“(D) satisfying 

A(—— Aijj) = in D] fj = ^ = 0 on dD, 

n — 1 av 


and 

(3.21) 


11^112+.. < Cplieil-., 


s = 0,1 
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where ri € (0,1], tq € (0, 2], Cp denotes the prior constant dependent on n{x) and D 
but independent of the right-hand side ^ of the equation. 

It is easy to see that (I3.2ip is valid with = 2—s when n and dD are appropriately 
smooth. When n is a constant and C is a convex polygon, from Theorem 2 in 
[3] we can get ri = 1 and if the inner angle at each critical boundary point is smaller 
than 126.283696...° then tq = 2. 

Consider the auxiliary boundary value problem: Find G Hq{D) such that 

(3.22) A( - Acj)i) = —A(u — Piu), in D, 

n — 1 

(3.23) </,! = ^ = 0, on dD. 

ou 

Let R(D) hold, then 

Ulh+r, < ||A(m - PhU)\\-i < ||W - 
Consider the auxiliary boundary value problem: 

(3.24) A{^—-A(j)2) = u-P},u, inD, 

n — 1 

(3.25) (j )2 = = 0, on dD. 

dn 

Let R(D) hold, then 

||</'2||2-hro < \\u-Plu\\o. 

Lemma 3.4. Suppose that R(D) and ( C2) are valid (s = 0,1), then for {u, to) € H, 

(3.26) \\{u,uj) - Phiu,uj)\\H‘' <h''‘\\{u,uj) - Phiu,uj)\\ii, s = 0,l. 


Proof The weak form of (I3.22|1 - (I3.23|1 is 

(A^;, ^A<^i)o = (V^;, V(m - PIu))o, G H^{D). 

Let V = u — P^u, then 

(3.27) (A(zr - P^u), = {V{u - P^u), V{u - P^u))o, 

which leads to 

l|V(u - Plu)\\l < ||A(w - P^)||o||A(0i - P^Vi)l|o 
< \\A{u - Plu)\\oh^^\\u - P^uWi, 

thus 

(3.28) \\u-Pl,uh<h''-\\u-P^u\\2. 

The weak form of (|3.24I) - (I3.25I) is 

(A?;,--^A^2)o = ('i^,M - VvGHq{D). 

















ERROR ESTIMATES FOR TRANSMISSION EIGENVALUES 


11 


Let V = u — P^u, then 

(3.29) (A(w - Plu), -^—A{cj)2 - Phh))o = {u - PIu,u- Phu)o, 

n — 1 

thus 

\\u-Phu\\l < \\A{u - Pj[u)\\o\\A{(j)2 - Ph(l)2)\\o < \\u - Pluhh’^'^Wu - P^uWo. 

and 


(3.30) \\u-P^u\\o<h-o\\u-P^u\\2. 
Combining (13.281) and (I3.30p we get 

(3.31) \\u-P^u\\s ^h^-’Wu-Plu\\2, s = 0,1. 
By the definition of the negative norm, 

(w-P^u;,Oo 


W - PhUj\\s -2 = 


(3.32) 


< 


sup 

^o~ 


sup 

n- 


llel| 2 -. 

u;-Pioj\\oU-Piao 


llel| 2 -. 


< h^-‘^\\oj - P?UJ 




s = 0,1. 


Ilw - 


h^Wo 


Combining (13.311) with (I3.32|) we deduce that 

||(M,w)-P/,(u,a;)||H« = \\{u - Plu,uj - P^u;)\\h‘> 

< ||n - Plu\U + llo. - Plu;\U_2 < h^^Wu - P^uh + h^-^\ 

< + /i^“'*)||(m, w) - P/i(m,w)||h, s = 0 ,1, 

noting that h?~'^ < we obtain (13.261) . □ 

Theorem 3.5. Suppose that R(D) and ((72) are valid, and n G W^~^’°^(D) 
fs = 0,1). Then 

(3.33) d{R{E),RiEh))s<h^^eh{X), s = 0,1. 

Let {uh,0Jh) with \\(uh,uJh)\\A = 1 be an eigenfunction eorresponding to Xh, then there 
exists eigenfunetion (u,ui) eorresponding to X, sueh that 

(3.34) ||(m?»,w/i) - (m, w)||h- < (^'''’e/i(A))“, s = 0,1; 
further let a = 1, then 

(3.35) \\{uh,uJh) - {u,uj)\\h‘> < h'^‘’\\{uh,w!h) - {u,uj)\\h, s = 0,1. 


Proof. Note that in R{E) the norm || • Hh® is equivalent to the norm || • ||h, and 
TR{E) C R{E), by (I3.26|) we deduce that 

\\T{u,u}) - Th{u,u})\\ii- 


||(P-7)i)|r(e)||h« = sup 

{u,Lj)eR{E) 




< 


||r(it,a;) - PhT{u,u)) 




sup 

{u,Lj)eR{E) IK'u,a;j||H 

\\T{u,uj) - PhT{u,uj)\\ii 


< /l’'® sup 

{u,uj)gR{E) 




(3.36) =V- sup l|r(».^)-Psn“.^)"H 
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Thanks to Theorem 3.1, Theorem 7.1 in [2] and (13.361) we deduce 
9{R{E),R{Ef,)),<\\{T-n)\R^E)\\n> <h^^en{X). 
And from Theorem 7.4 in we get 

\\{uh,uJh) - (w,w)||h« < ll('r-r,,)|/j(£)||g, < (/i’'^e/,(A))°. 
When a = 1, we have 


\\{Uh:UJh) - {u,Uj) 


H' 


< 


||(r-T^)(u,a;) 




< 


(u,a;) - Ph{u,uj) 




thus from (I3.26|) and IKwjw) — Ph{u,uj)\\H. ^ ll('w,w) — {uh,ojh)\\H we obtain (13.351) . 
This completes the proof. □ 

Remark 3.1. Using the same argument as in this section we can prove the error 
estimates of finite element approximation for the dual problem (12.151) . for example, 
there hold the following two estimates 


(3.37) 

(3.38) 


«,<)-(«*,cu*)||h<4(A*)^ 


\ / * * 
- (m 




< 


s = 0,l. 


4. T-wo grid discretization scheme. In this section we use the two grid dis¬ 
cretization scheme to treat transmission eigenvalues problem. 

Definition 4.1. V (u, z), (v*, z*) G H, B{{v, z), (v*, z*)) ^ 0, define 
as the generalized Rayleigh quotient of {v,z) and {v*,z*). 

We now outline the two grid discretization scheme. 

Scheme 4.1. Two grid discretization scheme. 

Step 1. Solve (13.11) on a coarse grid tth- Find Xh G C, {uh,ujh) G H// such that 
\\{uH,iOH)\\A = 1 and 


(4.1) A{{uh,ujh),{v,z)) = XhB{{uh,ujh),{v,z)), V(i;,z)gH//, 


and find (u^,wjj) G Mh{X*) according to Lemma 4.1 and Remark 4.1. 

Step 2. Solve two linear boundary value problems on a fine grid ■nh'- Find 
G H/i such that 

A{{u^,uj'^),{v,z)) = XhB{{uh,ujh),{v,z)), V(r;,z) G H^, 
and find G H/i such that 

A{{v,z),{u^*,u'^*j) = XHB{{v,z),{u*H,0J*H)), y{v,z)GHh. 

Step 3. Compute the generalized Rayleigh quotient 


A basic condition to the two grid discretization scheme is that \B{{uh(uh-i ^|f))l 
has a positive lower bound uniformly with respect to 77 (see, e.g.. Theorem 3.5 in [24]). 
The following Lemma 4.1 guarantees this condition. Using this condition, in the fol¬ 
lowing Theorem 4.3 we can also prove \B{{u^,u!^*))\ has a positive lower 
bound uniformly with respect to h. 

Lemma 4.1. Suppose that Xh = Xj^h (j = k,k + 1, ■ ■ ■ + q — i), and {uh,ojh) 

is an eigenfunction corresponding to Xh- Let he the orthogonal projection 

of {uh,<jJh) to Mh{X*) in the sense of inner product A(-, and let 

{Uh, OJh) 

\\{uh,^h)^a 


(4.2) 


{u*H, UJ^) 
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Then when H is small enough has a positive lower bound 

uniformly with respect to H. 

Proof Define f{{v,z)) = A{E{v,z), (u',uj')), where = \\e(uh,’^h)\\a ' 

for any (v^ z) G H 

\f{{v,z))\<\\E[v,z)\\A<\\EU\\{v,z)U 

/ is a linear and bounded functional on H and ||/||a < ||£'IU- Using Riesz Theorem, 
we know there exists (u, w) € H satisfying 

^((u,z),(u,w)) = A{E{v,z),{u',uj')), ||(u,w)|U < ||£^IU- 

Now, for any (u, z) G H, 

A((v, z), (A*-i - TT(^,^)) = Al((A-i - T)“(u,z), (u,^)) 

= A(E(A-^ - r)“(u, z), (u', cv')) = A((A-^ - TrE(v, z), (u', lu ')) = 0, 

i.e., (A*-i-T*)“(iI,S) = 0. Hence (u,di) S R{E*). Let : 

then {u'e,uJe) £ R(Ee) = Mh{A*), 

1 = H((uh,wh), {u'h,uj'h)) 

= A{{uh,UJh) - {Uh,UJh) + {u'h.uj'h)) 

( 4 - 3 ) = A{{uJi,uje), {u'h,uj'h)) < \\{uh,uje)\\a\\{u'h,ujh)\\a- 

From [a [IQ], if iL ^ 0 we have 

\\{Ee - E*)\ji(E‘‘)\\A -A 0, 

\\E{uH,tVH)\\A = \\E{uh,UJh) — {uH,tVH) + {uh,ujh)\\a -a 1, 

and 


A{{uh,ujh),Ee{u,uj)) = A{{uH,iOH), {Efj - E*){u,uj)) + A{{uh,ujh),E*{u,uj)) 
= A{{uh,ujh), {E*e — E*){u,lo)) + \\E{uht^h)\\a -a 1. 


Thus, there is a constant Cq independent of El such that 


(4.4) 


\\{u'eM\\a = 


\E 


hk 


\A{{uh,ioh), Efj(u,uj))\ 


<Co- 


From (13.IL (14.31) and (14.41) . deduce 


\B{[uh,ujh),{uh^^h))\ 
= \^~h\\M{uh,^h), 






Uh^^h)\\a 

1 


■)l — 


|A 


> 


H\ 




A UolAi/l 


{Uh, uj^) 
\\{'^HtA>h)\\a 


Then when El is small enough \B{(uht<^h),{u*e,uj*e))\ has a positive lower bound 
uniformly with respect to H. U 

The above lemma is fundamental for studying two grid method as well as multigrid 
method and adaptive algorithm for the transmission eigenvalue problem and many 
2 mth order nonselfadjoint elliptic eigenvalue problems. 
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Remark 4.1. Computational method for 

Step 1. Find a basis {(V'A Mh{X*)- 

How to seek this basis efficiently is an important issue of linear algebra. When 
the ascent of Xi,H is equal to 1 {i = k, - ■ ■ ^k + q—1), we can use the Arnoldi algorithm 
[101 [miss] to solve the dual problem of 63) and obtain this basis and meanwhile 
MATLAB has provided implemented Arnoldi solvers “sptarn” and “eigs”; we can also 
use the two sided Arnoldi algorithm in m to compute both left and right eigenvectors 
of (14.11) at the same time, and obtain a basis {(uy//, in Mh{X) and a basis 

Step 2. Solve the following equations: find /3i€C, i = + 1 such that 

k+q-l 

X! l = k,---,k + q-l. 

i—k 

Then and = (u^, w^)/||(u^,a;;^)|U- 

How to determine the multiplicity g of A is a difficult task in mathematics and has 
few achievements by now. In practical computation, q is determined by computational 
experience. When q cannot be determined exactly, we also use q\g to replace g, where 
qxj^ is the multiplicity of Xh that can be easily determined in practical computation. 

Our analysis makes use of the following lemma which is a generalization of Lemma 
9.1 in |5] and Lemma 3.6 in [24] . 

Lemma 4.2. Let (X,u,uj) and (A*,m*,u;*) be the eigenpair of i2.6\} and (2.151) . re¬ 
spectively. Then, \/{v, z), {v*, z*) € H, B{(v, z), {v*, z*)) ^ 0, the generalized Rayleigh 
quotient satisfies 

A{{v,z), (v*,z*)) _ ^ ^ A{{v,z) - {u,uj), {v*,z*) - {u*,uj*)) 
B{(v,z),(v*,z*)) B{{v,z),{v*,z*)) 

rx _ . B{{v,z) - {u,uj),{v*,z*) - iu*,uj*)) 

^ ’ B{{v,z),{v*,z*)) 

Proof. From (12.61) and (12.151) we have 

A{{v,z) - {u,uj),iv*,z*) - -XB{{v,z) - (u,oj), (v*, z*} - (u*,aj*)) 

= A((v,z),(v*,z*))-A((v,z),(u*,uj*)) 

- A((u,uj), (v*,z*)) + A((u,uj), (u*,uj*)) 

- MB((v,z), (v*,z*)) - B((v,z), (u*,uj*)) 

- B((u,uj), (v*,z*)) + B{{u,uj), {u*,uj*))) 

= A{{v,z), {v*,z*)) - XB{{v,z), iu*,uj*)) 

- XB{{u,uj), iv*,z*)) + XB{{u,uj), {u*,uj*)) 

-X{B{{v,z),{v*,z*))-B{{v,z),{u\lv*)) 

- {v*,z*)) + B{{u,uj), {u*,uj*))) 

= Aiiv, z), {v*,z*)) - XB{{v, z), iv*,z*)), 

dividing both sides by B{{v, z), {v*,z*)) we obtain the desired conclusion. □ 

Theorem 4.3. Let Xh, {uh, ojh), (u*ht'^*h)^ be the numer¬ 

ical eigenpairs obtained by Scheme 4-1, o-n-d Xh = Xj^h (j = k,k + 1, ■ ■ ■ ,k + q — 1). 
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Assume that the ascent of A are equal to 1, and that R(D) and {C2) are valid, 
n S (s = 0, IJ. Then there exists {u,uj) G R{E) and {u*, uj*) G R{E*) 

such that when El is properly small there hold 

(4.6) \\{u\lu>^) - («,w)||h < H-^enW+ehiX), 

(4.7) ||(u"*,a;"*) - («*,w*)||h < + el{X), 

(4.8) I A'‘-A|<{i7'-=ejf(A)+e,,(A)}{iJ’-=£j,(A*)+4(A*)}, s = 0,l. 


Proof. Let (u,a;) G R{E) such that {u,uj) — {uh,(^h) and A — Xh both satisfy 
Theorem 3.3 and Theorem 3.5. From 1)2.121) we get {u,uj) = XT{u,uj), and from (13.31) 
and Step 2 of Scheme 4.1, we get {u^,uj^) = XHTh{uH,aJH)- From p.3l) . we have 

\\XHTh{uH,ajH) — Ar;^(u,a;)||H < ||A_y(u//,w//) — A(m,w)||h'>, 

and thus by using Theorem 3.3 and Theorem 3.5, we derive that 

||(u^,a;^) — (u,a;)||H = \\XHTh{uH,a:H) — AT(u,a;)||H 

< ||AffT?i(u_y,w_y) — XTh{u,Lo)\\H + \\XTh{u,uj) — AT(u,a;)||H 
^ ||Aff(u_f/,w_y) — A(u,a;)||H'> + |A|||I)i(M,w) — T(u, w)||h 
<i7"“eff(A)+e;,(A), 

i.e., (14.61) holds. Similarly we can prove (14.71) . 

From (14.5L we have 

. _ A{{u\ u;^) - (u, u;), - {u*,u;*)) 

(u^*,w^*)) 

i3((u^cc'-) - (u,cc), (u^*,cc^-) - {U*,UJ*)) 

Note that {uhtOJh) and (it^,w^) just approximate the same eigenfunction (u,w), 
and (u^*,w^*) approximate the same adjoint eigenfunction {u*,u!*), and 
\B{{uH,ajH), {u*ij,u!'^))\ has a positive lower bound uniformly with respect to H. 
From 

i?((u^a;'‘), (u"*,a;^*)) = B{{u\uj'^), - B{{u,u), iu*,u;*)) 

+ B{{u,u}), iu*,uj*)) - B{{uh,ujh), {uh,ujh)) + B{{uh,ujh), {uh,wjh)), 

we know that |i?((M^, w^), {u^*, uj^*))\ has a positive lower bound uniformly. Therefore 
from (14.9L we get 

(4.10) lA'^ - A| < ||(u^w'‘) - {u,co)\\u\\{u>^*,uj>^*) - (w*,w*)||h. 

Substituting (14.611 and (14.711 into (14.1011 yields (14.8F □ 

5. Numerical experiment. In this section, we will report some numerical ex¬ 
periments for the finite element discretization m and two grid discretization scheme 
(Scheme 4.1) to validate our theoretical results. 

We use MATLAB 2012a to solve (I1.1F()1.4I1 on a Lenovo G480 PC with 4G mem¬ 
ory. Our program is implemented using the package iFEM m- 
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Let be a basis of and Uh — Ui^i,uJh = Stfi ^iii- A similar defini¬ 
tion can be made for and ujI^. Denote it = (ui, • • •, and at = (a;i, • • •, 

Similarly wt and wt can be defined from and uj'^. To describe our algorithm, we 
specify the following matrices in the discrete case. 


Matrix 

Dimension 

Definition 

Ah 

Nh X Nh 

“ Jd n—l 

Bh 

Nh X Nh 

Ki = + • niAdx 

Ch 

Nh X Nh 

— ~ Jd n — l 

Dh 

NhXNh 

fr> 


when n € W^’^{D), 6,,, = /^{V(;^6) ' 

Then (13.11) and (13.111) can be written as the generalized eigenvalue problems 


(5.1) 


(Ah 0 \ ^ 

0 Dh J\ it 




Bh Ch 
Dh 0 


it 

it 


and 

(5.2) 


(Ah 0\(u*\_ (BhDh 
0 Dh J J -^'^[Ch 0 


U* 

ZZJ* 


Note that in (EH), Ah is a positive definite Hermitian matrix and Dh can be 
equivalently replaced by the identity matrix Ih, which will lead to two sparser coef¬ 
ficient matrices with a good properties. Thus the computation of the eigenpairs for 
(EH) is efficient. Similarly, the matrices for Scheme 4.1 can be given but are omitted 
in this paper. 

For convenience, we use the following notations in our tables and figures: 

The jth eigenvalue obtained by EU) on TTh- 

kj,H = The jth eigenvalue obtained by (13.11) on -kh- 

k( = The jth eigenvalue obtained by Scheme 4.1. 

—: Failure of computation due to running out of memory. 


5.1. Model problem on the unit square. We first consider the case when D 
is the unit square [0,1]^ and the index of refraction n = 16 or n = 8 -I- cci — a: 2 . We use 
BFS element to compute the problem, and the numerical results are shown in Tables 
5.1-5.2. We also depict the error curves for the numerical eigenvalues (see Figure 5.1). 

According to regularity theory, we know rg = 2 and u,ui G H'^{D). When the 
ascent a = 1: according to (13.191) and (13.151) . the convergence order of the eigenvalue 
approximation kj^h is four; according to (14.8|) . when h > H^, the convergence order 
of the kj is also four, i.e.. 


(5.3) 


\kj-hh\<h\ \k,-k^\<h\ 


It is seen from Figure 5.1 that the convergence order of the numerical eigenvalues 
on the unit square is four, which coincides with the theoretical result (15.31) . 

In addition. Tables 5.1-5.2 show that the numerical eigenvalues obtained by BFS 
elements give a good approximation; it is worthwhile noticing that two grid discretiza¬ 
tion scheme can achieve the same convergence order as solving the eigenvalue problem 
by BFS element directly. Moreover, we find that the two grid discretization scheme 
can be performed on finer meshes so that we can obtain more accurate numerical 
eigenvalues. Therefore, the two grid discretization scheme for solving transmission 
eigenvalue problem is more efficient. 
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Table 5.1 

The eigenvalues obtained by BPS element on the unit square^ n = 16. 


3 

H 

h 


k- 

kj,h 

1 

x/2 

8 

V2 

32 

1.8800518272 

1.8795932933 

1.8795931085 

1 

16 

V2 

128 

1.8796216444 

1.8795911813 

1.8795911812 

1 

V2 

32 

V2 

256 

1.8795931085 

1.8795911697 

1.8795911747 

2 

V2 

8 

v/2 

32 

2.4462555154 

2.4442475976 

2.4442447101 

2 


v/2 

128 

2.4443713201 

2.4442361446 

2.4442361333 

2 

32 

x/2 

256 

2.4442447101 

2.4442361002 

— 

3 

8 

t/2 

32 

2.4462555154 

2.4442475976 

2.4442447101 

3 

v/2 

16 

n/2 

128 

2.4443713201 

2.4442361446 

2.4442361334 

3 

V2 

32 

n/2 

256 

2.4442447101 

2.4442361002 

— 

4 

8 

32 

2.8681931483 

2.8664515120 

2.8664469634 

4 

22 

V2 

128 

2.8665606968 

2.8664391607 

2.8664391408 

4 

256 

2.8664469634 

2.8664391111 

— 


Table 5.2 

The eigenvalues obtained by BPS element on the unit square, n = 8 + xi — X 2 , 


3 

H 

h 


k- 

kj,h 

1 

V2 

8 

y/2 

32 

2.8234457937 

2.8221946996 

2.8221945051 

1 

V2 

16- 

V2 

128 

2.8222709846 

2.8221893629 

2.8221893619 

1 

32 

vi 

256 

2.8221945051 

2.8221893480 

— 

2 

v/2 

8 

V2 

32 

3.5424522436 

3.5387180040 

3.5387126105 

2 

V2 

16. 

V2 

128 

3.5389469722 

3.5386967795 

3.5386967579 

2 

32 

\/2 

256 

3.5387126105 

3.5386967008 

— 

5,6 

8 

32 

4.4971031374 

±.8770188489i 

4.4964665266 

±0.87155340261 

4.4965591247 

±0.8715053132i 

5,6 

32 

V2 

128 

4.4965591247 

±0.8715053132i 

4.4965519816 

±0.8714818735i 

4.4965519832 

±0.8714818728i 

5,6 

32 

256 

4.4965591247 

±0.8715053132i 

4.4965519531 

±0.871481788i 

z 


5.2. Model problem on the L-shaped domain. We consider the case when 
D={—1, 1)^\([0,1) X (—1, 0]) is the L-shaped domain and the index of refraction n = 16 
OT n = 8 + Xi — X 2 - The numerical results obtained by BPS element are shown in 
Tables 5.3-5.4 and Figure 5.2. 

Figure 5.2 indicates that on the L-shaped domain, the convergence order of 
ki,h,k 2 ,h is around 1.3 and 2.3 respectively for n = 16, and the convergence order 
of ki^h,k 5 ^h is around 1.4 for n = 8 -f xi — X 2 - This fact suggests that the four 
eigenfunctions on the non-convex domain do have singularities to different degrees. 

5.3. Model problem on the circular domain. We also investigate the case 
of n being piecewise constant for a disk D of the radius 1. Let n(x) = ni for x G Di 
and n(x) = n 2 for x G D\Di, Di being an inner disk of the radius ri < 1. For 
this disk domain we generate a triangular mesh with h « ^ and number of degrees 
of freedom 101040 and move the nodes outside and nearest the inner disk onto the 
inner circle. And we use the Argyris element to solve (I1.1D - (I1.4I) on the mesh and the 
numerical eigenvalues associated with different n and ri are shown in Table 5.5. The 
analytic eigenvalues can be referred to Tables 2 and 3 in [19]. 
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Table 5.3 

The eigenvalues obtained by BFS element on the L-shaped domain, n = 16. 


3 

H 

h 


k- 

kj,h 

1 

8 

32 

1.4850653844 

1.4781249432 

1.4780403370 

1 

%/2 

16 

64 

1.4802422297 

1.4770298927 

1.4770116105 

1 

vT 

128 

1.4802422297 

1.4765529659 

— 

2 

V2 

8 

n/2 

32 

1.5705634174 

1.5697720130 

1.5697716222 

2 

8 

n/2 

64 

1.5699010557 

1.5697385580 

1.5697385335 

2 

V2 

128 

1.5699010557 

1.5697293878 

— 

3 

32 

1.7078128918 

1.7055170229 

1.7055794458 

3 

V2 

16 

V2 

64 

1.7061981331 

1.7052820968 

1.7052949545 

3 

v/2 

16 

n/2 

128 

1.7061981331 

1.7051443044 

— 

4 

8 

vT 

32 

1.7834680505 

1.7831208381 

1.7831208523 

4 

%/2 

m 

V2 

64 

1.7831489498 

1.7831171004 

1.7831171097 

4 

n/2 

128 

1.7831489498 

1.7831163182 

— 


Table 5.4 

The eigenvalues obtained by BFS element on the L-shaped domain, n = 8 xi — x^- 


3 

H 

h 


fc? 

kj,h 

1 

\/2 

8 

V2 

32 

2.3127184233 

2.3045393229 

2.3043808707 

1 

vT 

16 - 

64 

2.3069612717 

2.3032153953 

2.3031811119 

1 

16 . 

32 

2.3069612717 

2.3026641395 

— 

2 

8 

2.3974892428 

2.3957871378 

2.3957863810 

2 

16 . 

V2 

64 

2.3960567236 

2.3957182671 

2.3957182148 

2 

16 - 

V2 

128 

2.3960567236 

2.3956994585 

— 

5,6 

8 

V2 

32 

2.9287346709 

2.9252416365 

2.9257099237 


±0.5743458101i 

±0.56645452481 

±0.56643078811 

5,6 

16 


2.9272401495 

2.9248351944 

2.9249318006 



±0.56860457021 

±0.56545378851 

±0.56544872001 

5,6 

16 

V2 

128 

2.9272401495 

2.9244335090 

— 


±0.56860457021 

±0.56499946781 

— 




Fig. 5.1. Error curves on the unit square for k\,k 2 with n = 16 (left), and for k\,k^ with 
n = 8 ^ xi — X 2 (right). 
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Fig. 5.2. Error curves on the L-shaped domain for ki, k 2 with n = 16 (left), and for ki, k^ 
with n = S xi — X 2 (right). 


Table 5.5 

The eigenvalues obtained by Argyris element on the circle of radius 1 for piecewise constant n. 


ni,n2,ri {ki.h’ ^2.h, k^.h’ ^6.h, h.h, kg^h, ho.h,/^ll,h} _ 

13 5 0.5 1.4975,1.7353,1.7357,2.1709,2.1721,2.3570±0.4881i,2.6995,2.6999 

2 4 0.5 2.4349±0.6969i,3.5812±0.5731i,3.5824±0.5730i,3.9137,3.9196,4.1160 

5 8 0.6 1.7891,2.2512,2.2520,2.5574±0.3780i,2.6681,2.6717,3.0376,3.0388 

5 3 0.6 2.2996±0.7437i,2.8876,2.8885,3.1861,3.2953,3.2982 

10 8 0.7 1.3732,1.7322,1.7328,2.1101,2.1121,2.3665±0.4401i,2.4958,2.4964 

2 4 0.7 2.4280±0.6299i,3.7582±0.6135i,3.7593±0.6136i,4.9689±0.4384i, 4.9751±0.4395i,5.0515 

13 11 0.8 1.1499,1.4899,1.4903,1.8261,1.8278,2.1570,2.1575 

3 6 0.8 2.2223±0.4796i,3.0286,3.0297,3.6711,3.6748,3.8012 

6 13 0.9 1.5929,2.0271,2.0280,2.4687±0.2643i,2.5142,2.5171,3.0034,3.0043 

6 2 0.9 2.0222,2.4190,2.4192,2.7401ili0.4468i,2.9309,2.9317 
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